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Abstract The covariant Klein-Gordon equation requires twice the bound- 
ary conditions of the Schrodinger equation and does not have an accepted 
single-particle interpretation. Instead of interpreting its solution as a proba- 
bility wave determined by an initial boundary condition, this paper considers 
the possibility that the solutions are determined by both an initial and a final 
boundary condition. By constructing an invariant joint probability distribu- 
tion from the size of the solution space, it is shown that the usual mea- 
surement probabilities can nearly be recovered in the non-relativistic limit, 
provided that neither boundary constrains the energy to a precision near h/to 
(where to is the time duration between the boundary conditions). Otherwise, 
deviations from standard quantum mechanics are predicted. 

Keywords Relativistic Quantum Mechanics • Classical Fields • Quantum 
Foundations 
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1 Introduction 

Because the Klein-Gordon equation (KGE) for a classical, scalar field was 
the conceptual precursor to Schrodinger's "derivation" of the Schrodinger 
Equation (SE), there is a widely- held misconception that the SE must be the 
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non-relativistic limit of the KGE. In fact, a second-order (in time) differen- 
tial equation like the KGE does not reduce to a single first-order differential 
equation (like the SE) in any limit One can reduce the KGE to a single SE 
only by artificially discarding half of the solutions - the so-called "negative- 
energy" solutions that evolve like exp(+icot). This may seem reasonable when 
considered in the framework of standard, non-relativistic quantum mechanics 
(NRQM), but NRQM must be a limit of a generally covariant theory where 
such a step is not natural. Indeed, this historical division between two arbi- 
trary halves of the solution has created severe problems for quantum gravity, 
for in curved space-time there is no well-defined way to separate the positive- 
and negative-energy terms pQ . 

This paper motivates and analyzes a novel interpretation of the KGE 
(corresponding to a neutral, spinless particle) that promises to reproduce 
the predictions of standard NRQM without arbitrarily discarding any of the 
solutions. Efforts to give the KGE a probabilistic interpretation analogous 
to the SE have historically led to two problems: apparent negative-energy 
solutions and no positive-definite probability four-current. Several recent pa- 
pers have addressed these issues [2] [3] [4] [5] , but there is not yet any accepted 
resolution. Both of these problems arise from the presumption that any in- 
terpretation of the KGE must be directly analogous to the standard interpre- 
tation of the SE. However, such a direct analogy is inherently unlikely, due 
to the mathematical disconnect between first-order and second-order differ- 
ential equations. Specifically, while the full solution to the SE requires only 
an initial boundary condition ijj(t = 0), the full solution to the KGE requires 
two independent initial boundary conditions: both <p(t = 0) and the first time 
derivative <p{t = 0). Without some expanded measurement theory to explain 
how to simultaneously impose independent boundary conditions on both <f> 
and 4>, one cannot even solve the KGE, let alone interpret the solutions. 

But if there was some way to impose additional initial boundary con- 
ditions (outside the scope of standard QM), then that extra information 
should allow for experimental predictions that surpass QM. A century of 
experimentation has not accomplished this feat, providing strong evidence 
that one cannot physically impose the necessary initial conditions to solve 
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the KGE. In other words, for a generic (neutral, spinless) particle, the maxi- 
mum amount of knowable information can be encoded by the instantaneous 
values of a single complex scalar field ip - only half as much information 
as the instantaneous values of the <fi and </> fields needed to solve the KGE 
for a complex scalar field. This "half-knowledge" situation is curiously simi- 
lar to the axiomatic foundation of Spekkens's interesting toy model [Sj, but 
still leaves the unanswered question: how can one calculate the solutions to 
a second-order differential equation like the KGE without sufficient initial 
boundary conditions? 

I propose that a natural resolution to this dilemma can be found by con- 
sidering time-symmetric approaches to NRQM; namely, those in which one 
imposes two boundary conditions at two different times corresponding to a 
preparation and a measurement 7J[8J[9J. If only half of the required infor- 
mation to solve the KGE can be imposed as an initial boundary condition, 
then a subsequent measurement can still impose the remaining constraints 
as a final boundary condition. This approach can also be extended to curved 
space-time by imposing a closed hypersurface boundary condition on clas- 
sical fields (similar to the closed boundaries imposed on quantum fields in 
recent work by Oeckl [TU] [H] ) . 

To compare this approach to NRQM, one simply constrains the KGE 
by imposing two boundary conditions on two different instantaneous hyper- 
surfaces (corresponding to consecutive external interactions/measurements). 
The solution to the KGE can then be "retrodicted" in the space-time volume 
between the two measurements. Because the full solution cannot be deter- 
mined by the initial measurement alone, one is forced into a probabilistic 
description, weighting the solution space to determine the relative likelihood 
of different pairs of measurements. After the final measurement, when the 
additional boundary conditions become known, such a probabilistic interpre- 
tation is no longer necessary. The central result of this paper is to demonstrate 
the existence of a relativistically invariant weighting scheme that, in the non- 
relativistic (NR) limit, gives probabilities very similar (but not identical) to 
NRQM. Furthermore, the historical problems with the KGE are solved by 
this approach: all possible solutions to the KGE represent positive energy 
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(as is already the case for a classical scalar field), and the probabilities are 
always positive definite. 

This paper is organized as follows: Section 2 solves the two-boundary 
Klein-Gordon equation for generic time-even boundary conditions, circum- 
venting the infinite poles in the same manner as in quantum field theory. Sec- 
tion 3 then motivates an invariant joint probability distribution, and demon- 
strates that it is nearly identical to standard conditional probabilities in the 
NR limit (given other constraints that are derived in an appendix) . Section 4 
addresses time-odd, multiple, and incomplete measurements. Section 5 exam- 
ines the constraints from the appendix, and finds that they are closely related 
to a time-energy uncertainty principle. Section 6 summarizes the postulates 
used in this approach, details some of the next steps required to extend this 
research program, and touches on other implications. 



2 Two-Time Boundary Conditions 

In the absence of a potential, the Klein-Gordon equation on a complex scalar 
field 4> (and its solutions) can be written 

cj,(r, t)= J a(fe)e i(fc -''^ t) + b(k)e i(k - r+UJt) dk. (2) 

Here a(fc) and b{k) are complex functions, and to is assumed to be a positive 
function of k; uj(k) = y c 2 k 2 + m 2 c' i /h 2 . 

Solutions of the form ([2]) have exactly twice the independent parameters 
as the solutions to the Schrodingcr Equation (for which the b(k) terms are all 
identically zero). The analysis from the previous section can equally well be 
applied to this doubled parameter space: If an initial preparation constraint 
can only specify parameters equivalent to a(k) (a single complex function), 
then specifying the two independent complex functions in the KGE solutions 
must require twice the constraints. Placing a second constraint at the next 
measurement is independently motivated by CPT-symmetry [9], but this 
symmetry implies both boundaries should be mathematically equivalent - 
not a strict initial boundary followed by some discontinuous "projection" at 
the final boundary. 
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There have been several prior efforts to impose two-time boundary con- 
ditions in NRQM without doubling the parameter space of the Schrodinger 
equation p~2] [13] ; these efforts have revealed that the equations generally be- 
come overconstrained. In light of this difficulty, other two-boundary efforts 
have used a doubled-parameter space [7,8,14], where the ordinary quantum 
wavefunction ip is constrained by an initial boundary and some other wave- 
function is constrained by a final boundary. In the present work, no artificial 
division is made between any two halves of the full Klein-Gordon field. (Re- 
call, it is this arbitrary splitting into so-called positive-energy and negative- 
energy solutions that causes trouble in curved space-time.) Instead, this ap- 
proach uses only one field, 0, all of which is constrained by both an initial 
and a final boundary. (Using this same approach, Rovelli has calculated the 
action of a classical scalar field [T3] , but no direct use of this result has been 
made except for a noted similarity to the quantum field propagator.) 

In order to impose boundary conditions, it is convenient to use conven- 
tional measurement theory from NRQM, where a measurement constrains <f> 
to be an eigenfunction of some operator at a particular time (on a space-like 
hypersurface s). This permits a boundary to be imposed on, say, <fi(r, t = 0), 
but only in the reference frame somehow defined by the measurement it- 
self. (This will eventually have to be generalized, as discussed in section 
6.) Another complication is that the momentum operator in physical space 
P — ► —ihV has been constructed for the SE, not the KGE. Applying it to 
solutions of the form {2j) leads to unphysical results, because the b(k) terms 
physically propagate in the opposite direction as determined by the eigenval- 
ues of P. This is a general problem for any time-odd operator applied to the 
KGE. A natural solution will be presented in section 4; for now the analysis 
is simply restricted to time-even measurement operators such as X and P 2 . 

Suppose the eigenfunction of the initial time-even measurement operator 
on the hypersurface t = is f(r), and the eigenfunction of the final time-even 
measurement operator on the hypersurface t = to is g(r). (These eigenfunc- 
tions are determined via measurements in precisely the same manner as in 
standard NRQM.) Fourier-expanding (and ignoring the 27r's, which will be 
automatically incorporated into a later normalization) leads to the general 
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boundary conditions 

0(r,O)=/(r) = J F(ky k - r dk, (3) 

<f>(r,t )=g(r) = J G{k)e lk ' r dk. (4) 

Comparing these equations to l[2]). the uniqueness of the Fourier transform 
implies 

a(fe) + b{k) = F{k), (5) 

a(fe)e~^ t0 + b{k)e luJt ° = G(k). (6) 

This is the correct number of equations to solve for the coefficients a(fe) 
and b(k), but there is a problem for the particular values of k at which 
cj(k n )to = nir. At these discrete values, solving ([5]) and ([6]) yield expressions 
for a(k) and b(k) with infinite poles (for arbitrary boundaries F and G). 

But quantum field theory is no stranger to infinite poles - the propagator 
for the KGE contains a term (E 2 — p 2 c 2 — m 2 c 4 ) -1 which blows up on the 
mass shell. The standard prescription for solving this problem is to make the 
assignment m 2 — > m 2 — it. After performing the integral, one takes the limit 
e — > 0. Because this is known to give acceptable results, it seems reasonable to 
use precisely the same method here. Implementing this change in m 2 changes 
the solutions to the KGE; 

poo 

c/)(r,t)= a(k)e i( - k - r ~^e- et + b(k)e l{k - r+wt) e" t dk. (7) 



where, as usual, all positive constants are absorbed into e itself. (Al/w term 
has also been absorbed into e, which is approximately constant in a NR 
approximation.) One concern raised by the form of J7|) is that it appears to 
diverge as t — > ±oo. But recall that we are imposing both an initial and a 
final boundary condition, so between these boundaries <fi remains finite (as 
we eventually will take e — * 0) . 

Now it is possible to solve the two boundary problem, using the new form 
of the solution (JT]). The original boundary conditions ([5]) and (J6]) now appear 
as 

a(k) + b(k) = F(k), (8) 
a(k)e- iuJto e- eto + b(k)e iwto e et ° = G(k). (9) 
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These are now exactly solvable for any two generic boundary conditions, 
yielding 

F(k)e tuJto e et ° - G(k) 
a{k) = ■ \ ' ■ , , , (10) 

= G(k)- F(fc)e-"'°e-'^ 

Notice that the addition of the small quantity e now prevents the poles at 
u>(k n )to — nir from diverging. 

Before proceeding, it is worth noting that this technique has already man- 
aged the feat of finding a continuous field </>(r, t) everywhere between any pair 
of time-even measurements. No "projection" is needed to get from one exact 
measurement to the next; one simply uses (ITU1) and (fTTj) to solve J7]) exactly. 
Still, this cannot be done until after the second boundary becomes known, 
so this is a retrodiction, not a prediction. Also, some of the coefficients still 
diverge as e — > 0, making it difficult to ascribe any "reality" to <\>. Various 
physical motivations for a non-zero e will be discussed in section 6. For now, 
one can take the already-common "agnostic" viewpoint, treating <p as simply 
a mathematical tool for calculating measurement probabilities. In the next 
section, it is shown that this can be accomplished even if e — > exactly. 



3 Measurement Probabilities 

Given a two-time-boundary framework, the use of conditional probability 
needs to be revisited. In standard NRQM, some initial measurement (Qp,F) 
at time t=0 determines the state/wavefunction = 0). (In this notation, 
Qf is the experimenter's choice of measurement, and F is the particular 
outcome.) Evolving -0 to the time to of the next measurement (Qg,G), one 
traditionally uses the conditional probability P(G\ip(to), Qg)- But this time- 
directed process simply does not work for the KGE because (Qf,F) is not 
sufficient to determine the solution cj). Furthermore, this paper speculates that 
(Qc G) determines 4>(to), so one could hardly have a conditional probability 
of G in turn determined by <j){to). 

In this framework, one needs a time-symmetric implementation of proba- 
bility, closer in spirit to scattering amplitudes in quantum field theory. Such 
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a symmetry can be naturally achieved by treating joint probability distri- 
butions (JPDs) as the fundamental quantities from which conditional prob- 
abilities can be derived. This approach assigns a joint probability to every 
possible initial/final measurement outcome pair Ft, Gj and the time between 
them At = to. If there was some JPD given by J(Fi, Gj, At) associated with 
each possible pair of outcomes (for a given pair of measurement operators Qf, 
Qg) one could generate any particular conditional probability P(G% |i<2, At) 
via the standard normalization procedure: J{F%, G\, At)/ J^j J{^2, Gj,At). 
For operators with continuous eigenvalues, the sum becomes an integral over 
all eigenstates of the chosen operator Qg- Standard zero-potential probabil- 
ities in NRQM can be recovered from the JPD 

(12) 

For any given F and to, the conditional probability of measuring any 
function G is given by this same expression, Pq(G\F, to) = Jq(F, G, to), if one 
constrains the magnitude of F and G such that J \F\ 2 dk — J \ G\ 2 dk = 1. 
(Note these constraints are simply chosen for familiarity, and should not be 
interpreted as a separate normalization; changing the "1" to any constant 
would have no effect, due to the normalization in the standard J — > P proce- 
dure described in the previous paragraph.) Of course, the above expression 
for Jo is not relativistically invariant. The challenge for any relativistic ex- 
tension of quantum mechanics is to devise an invariant JPD that is equal to 
Jo in the appropriate limit. 

Note that JPDs are not probability densities, so the conditional proba- 
bilities derived from them need not satisfy a continuity equation. Indeed, to 
talk about the probabilities of a measurement at a time between t = and 
t = to is a contradiction in terms, because the next measurement happens at 
t = to by definition. A similar conclusion was reached in Oeckl's recent work 
on "general boundary" quantum field theory [lOj. 

In relativistic quantum mechanics, the charge density p(r,t) cannot rep- 
resent probability because it can be negative. Still, because of the usefulness 
of this density in relativistic quantum mechanics, this is a reasonable place 
to begin searching for an appropriate JPD. For the Klein-Gordon equation, 



Jo{F, G, to) 



J F(k)G*{k)e- tuJto dk 
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the charge density is given by 



P(r,t) 



2mc 2 




90 
Hi 



~~di 



) 



(13) 



and is real, although not invariant (it is one component of a four-current). It 
is not possible to make p invariant without some reference unit four-vector 
that defines the time direction. 

Fortunately, the boundary conditions define two hypersurfaces, S\ and 
S2- These surfaces have associated inward-pointing normal four- vectors r\\ 
and r)2 (pointing forward in time from the initial boundary and backward in 
time from the final boundary). In the limit comparable to NRQM, s\ and S2 
are the hyperplanes t = and t = to , so it is plausible that rj is the needed 
reference unit vector, but of course it is only defined on the hypersurfaces. 
Because of this constraint, p is not well-defined in the volume between the 
boundaries - but one can still integrate p along the hypersurfaces themselves. 
So, generalizing to any closed hypersurface defined by the two boundary 
conditions (e.g. the two infinite planes t = and t = to), the hypersurface 
integral 



is a scalar that is plausibly related to probability, motivated both by in- 
variance and known results from relativistic quantum mechanics. (In this 
notation, rf is the inward-pointing four-vector unit normal to the integra- 
tion surface and 9 M is the four-gradient with <9 = d/d(ct). Summation over 
the index fi is implied.) 

In the special case that the closed hypersurface is defined by the planes 
t = and t = to, W can be evaluated using (0. (Expand both <fi and <fi* as 
integrals in k and k', and then the spatial integral yields 2Tr8(k~ k'), leaving 
only a single integral in fc.) Dropping the overall constants, the simplified 
result is 




(14) 



W = J uj\a\ 2 (l - e- 2eto ) - cj\b\ 2 (1 - e 2et0 ) - Im [2eab* (l - e" 21 "* )] dk. 

(15) 
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Next, a(k) and b(k) can be written in terms of the boundary conditions F(k) 
and G(k) using (fT0| and ITTj) . To lowest surviving orders in e, this yields 

„ 7 /■ wei [|^| 2 + |G| 2 ] - 2u>et Re[FG*]cos(cjt ) + e«n(w*o)(...) „ 

W = / b 5-5 tt, r afc. 

(16) 

This does not vanish as e — ► because one gets a periodic delta function 
from the identity 



<^o \ -> 

S) ( 2 ti + sin 2 (ujt a ) 

u v ' n— 



lim 2 — . 9 ^ - : = VV(5(a;£o - n7r). (17) 



The sin(wio) term in the numerator of (|16p was not fully shown because it 
is odd over the delta function, and vanishes; the cos(ioto) term is either 1 or 
— 1 depending on n. Expanding dk in spherical fc-coordinatcs (k = \k\), one 
can use the relationship k dk = uj cLu to reduce the radial integral to a sum, 
leaving only a 2-D integral; 

f °° —~\ n uj 2 k 

w = / £ — d F »i 2 + i g »i 2 - 2i?e [ F « G «0 ( l8 ) 

^ n=n 

In this notation F n — F(k n ), G n = G(k n ), and k n = \k n \, where k n is the 
solution to u> n (k n ) = nir/to, and uq is the smallest value of n for which there 
is such a solution. 

I now propose that an invariant joint probability distribution can be con- 
structed from W according to 

J ={W max -W min f. (19) 

Here W max is the maximum value W can attain when varying any uncon- 
strained parameters; similarly W m in is the minimum value. Such a range 
parameter is common in probability theory: one way to interpret (I19[) is that 
any pair of boundary conditions permits a two-dimensional space of solutions, 
with the allowed range in each dimension determined by allowed values of 
W. Picking a solution at random, particular pairs of boundaries that have 
a larger solution space would be more likely than other pairs with a smaller 
solution space. 

Still, W has no range at all if the parameters F, G, and to are all exactly 
specified. But regardless of the precision of the measurements, standard mea- 
surement theory always allows for one free parameter: the unknown relative 
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phase between the two eigenfunctions /(r) and g(r). This relative overall 
phase, 9, does not affect the first two terms in (fT8|) . so for complete measure- 
ments these terms will not contribute to any calculation of J using ()19D . 

J must match the standard JPD in the NR limit (| 1 2|) , but Jo is given by a 
3-D integral over k-space, not a 2-D integral and a discrete sum. Still, the sum 
in (|18p approximates an integral if neither F(k) nor G(k) change very rapidly 
when the magnitude |fc| changes by k n+ i — k n ; the precise conditions are 
derived in the appendix and discussed in section 5. Given these constraints, 
the appendix demonstrates that the third term in (fl8|) can be excellently 
approximated as 



To examine whether or not J « Jo, one must now use the fact that 
both f(r) and g(r) were originally eigenfunctions of time-even operators 
(time-odd operators are considered in the next section). The eigenfunctions 
of non-degenerate time-even operators can be constrained to be entirely real 
if their relative phase exp(iQ) is added explicitly [16] . (For the degenerate 
case, there is no relative phase information, so one is still free to choose 
real combinations of the real eigenfunctions.) So, forcing f(r) and g(r) to 
be real by adding the relative overall phase exp(i9), the approximation (j2"0|) 
simplifies the JPD (JTHJ) into the expression 



This assumes that one is in the NR limit such that the in (fT5|) is roughly 
constant and can be pulled out of the integral (overall constants are irrelevant 
to the unnormalizcd JPD). Given the complete freedom of the unmeasurable 
quantity 9, (HJ) implies W max = - W min , so the range W max - W min is just 
2W max . 

Additional simplifications result from f(r) and g(r) being real; both 
Im[F(k)G*(k))cos(ujto) and Im[F(k)G* (k)]sin(ujto) are odd functions in k, 
so their integrals cancel out of both Jo (f!2| and J (|2Tj) . Extracting the It- 
independent portion ujq = mc 2 /h from ui{k) = ujq + ui%(k) further simplifies 




(20) 




(21) 
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these expressions to 



J = (^J Re[FG*]cos(uj 1 t a ) dkj + (J Re[FG*]sin(oj 1 t ) dkj , (22) 
J = Re[FG*]cos8 [cos(ujoto)cos(ujitQ) ~ sin(ujoto)sm(uJito)] dk 
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(23) 

While similar, these two expressions are unfortunately not equivalent if to 
is constrained to be a precise value. But is practically impossible to constrain 
<o to an accuracy less than uj^ 1 (~10 -20 seconds for an electron, and even 
less for more massive particles). No current experiments can determine the 
measurement time to a sub-attosecond accuracy, so W must be maximized 
over the completely unknown angle uioto as well as over 9. The latter maxi- 
mization is trivial, as one simply forces cos 8 = 1 in (f23|) . And with complete 
freedom of the angle woto? the maximum value of (AsinuJoto ± B coswoto) 2 
is just A 2 + B 2 . Using this fact, one can compare (|22|) and (|23|) to confirm 
that J a; Jo, to within the accuracy of the previous approximations. 

Note that this procedure only works if the precise value of to becomes 
part of the solution parameter space, and is not merely an unknown random 
quantity. In other words, if to is constrained to be some particular value - 
even if that value is completely unknown - then J is no longer necessarily 
equal to Jo. 

From a foundational perspective, using J is preferable to Jo , because J is 
relativistically invariant, unlike (|12[) . But theoretical preference is no longer 
the only issue; even given realistic uncertainties in to, the above approxima- 
tions ensure that J and Jo are not exactly the same, opening the door to 
an experimental differentiation between this approach and NRQM. J and Jo 
begin to diverge in the relativistic regime, and also when the approximation 
(|2H)) begins to fail as discussed in section 5. 

The other regime in which new effects would be expected is if the value 
of to could be experimentally constrained to a value much less than (jJq ■ 
Although such capabilities currently seem far out of reach, it is worth con- 
sidering a future experiment that might constrain tooto to be approximately 
a multiple of tt. In that case, the sm(woto) term in (|23|) can be ignored, and 
one finds that J reduces to only the first term in (|22|) . This term is the square 
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the real portion of the standard NRQM transition amplitude, and of course 
there are many deviations from NRQM that would occur if the imaginary 
part of the amplitude was ignored. 



4 Time-Odd, Multiple, and Incomplete Measurements 

The previous sections have assumed two consecutive measurements of time- 
even quantities, such as position, P 2 , or L 2 (angular momentum). Measure- 
ments of time-odd quantities leads to difficulties, because (as noted in Section 
2) the b(k) terms in @ physically propagate in the opposite direction as de- 
termined by the eigenvalues of P. Another problem, as noted in Section 1, is 
that there are two independent (but equally important) mathematical objects 
on which initial values can be specified: both <j> and the first time derivative 
<j>. Initial boundaries have historically only been imposed on the field itself 
because of the reliance on first-order differential equations, but when using 
the KGE one should consider the possibility of constraining the value of <j) 
as well as </>. 

Both of the issues in the previous paragraph are linked to a simple fact: the 
measurement of a time-even quantity must yield the same value under time- 
reversal of the entire system, while the measurement of a time-odd quantity 
(like momentum) must change sign under time-reversal. In spinless NRQM, 
time-reversal is accomplished via a complex conjugation, which is why (in 
the position basis) time-even operators are real and time-odd operators are 
imaginary. But complex conjugation does not time-reverse the KGE, only the 
SE, so using "i's" to distinguish time-even and time-odd operators must now 
fail. What does get a sign-change upon time reversal is 0, the very object 
that standard measurement theory does not address. It is therefore tempting 
to replace the "i n in time-odd measurement operators with —d/dt. While 
the operator units are now wrong, this is easily fixed by multiplying these 
time-odd operators by the natural unit of time lo^ 1 = h/(mc 2 ). This section 
will demonstrate that such a replacement solves the momentum eigenvalue 
problem while also allowing measurements to constrain <j) instead of <p. 
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Under the substitution i — ► — ui d/dt in all time-odd operators, the 
corresponding eigenvalue equation Q4> = q<fi then changes to become 

iQ4> = wo# (24) 

for any standard time-odd operator Q (and eigenvalue q corresponding to 
eigenfunction </>). 

In the limit where one expects the ordinary NRQM operators to be valid, 
the frequency u>{k) at the relevant values of fc are all approximately the 
same; lu ujq = mc 2 /%. Using this approximation, the time-derivative of ([2]) 
is simply 




-icj / a(fc)e l(fe ' r ^ t) - b{k)e i{k - r+UJt) dk. (25) 



Inserting this into (|24[) , one finds that the ordinary eigenvalue equation Q<j) — 
q<p is recovered when b(k) = (the standard NRQM condition). On the other 
hand, if a(k) = 0, the equation now has the opposite eigenvalue, — q. This is 
exactly the sign-reversal needed to solve the momentum eigenvalue problem 
noted above; if Q is —ihV, and yields a measurement eigenvalue hko at t = 0, 
then ([24| implies a boundary condition a(k) + b(— k) = 5(k — ko). Regardless 
of the relative phase or weight between a(fco) & n d b(—ko), the combination 
a(ka)exp(iko ■ r — itot) + b(—ko)exp(—iko ■ r + itot) physically propagates in 
a direction aligned with fco- (This would not have been the case if one was 
using the ordinary eigenvalue equation; the b(ko) term has a phase velocity 
opposite that of the b(— fco) term.) 

In general, (|24p implies that the measured real eigenfunction (in k-space) 
F(k) of a time-odd operator from standard NRQM should be imposed on the 
combination a{k)exp{~iujt) + b{—k)exp{iujt) . The section 3 conclusion that 
J ~ Jq therefore remains valid for two consecutive time-odd measurements; 
one need only change b(k) — > b(~k) starting at ([5]). 

It remains to show that a time-odd measurement followed by a time- 
even measurement (or vice-versa) also gives results consistent with NRQM. 
Suppose a time-odd measurement at t = yields a k-space measurement with 
(real) eigenfunction F(k), and a time-even measurement at t = to yields a 
measurement with (real) eigenfunction g(r), both within an arbitrary overall 
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phase. Then the new boundary conditions can be written in terms of F(k) 
and G(k) (the latter is calculated from g(r) according to ([4J) as 

a{k) + b{-k) = F(k), (26) 

a{k)e- iuJto e- eto + b{k)e lujt °e et ° = G(fc). (27) 

These equations have a solution different from (fT0|) and (fTT|) . Substantial but 
straightforward algebra yields 

F ( k ) e2iuJt ° e2et ° + F (~ k ) - G(k)e- iuJt °- et ° - G{-k)e iuJto+et ° 

a \ K ) ~ e 2iLut +2et _ e -2iujt -2et ' ^ ' 

G(k)e luJto+et ° + G{-k)e- lbJt °- et ° ~ F{k) - F(-k)e- 2iut °- 2et ° 

(29) 

Inserting these new values of a(k) and b(k) into (fT5|) yields a different value 
of W than given in (fTBjl . The only relevant terms are those with both F 
and G, because all of the terms like G(fe)G*(— fe) are independent of the 
relative overall phase 9 and do not contribute to variation in W. Putting in 
the relative phase explicitly (FG* — > exp(i6)FG*) one can simplify using 
the fact that for the real function in k-space F*(k) = F(k) and for the real 
function in position space G*(k) = G(—k). To largest surviving orders of e, 
this yields 

r ~uoet Re [F{k)G*{k) ( e 3iuJto+ie + e - luJto+l6 + 2e- luJt °- lS )] 
~ J eHl + sin 2 (2uoto) 

(30) 

As e — > 0, a periodic delta function arises according to (jTTJ) , but with 
twice as many poles as before: 8(ujt$ — mr/2). At each pole, exp(3iujto) = 
exp(—iujt ), but the value of this quantity picks up a factor of i at each 
consecutive pole, making both the Re(FG*) terms and the Im(FG*) terms 
important. Applying the same approximation (|20[) to both of these terms, 
one finds that in the NR limit 

W^cosO J Re[F(k)G*(k)]cos(ujto) ~ Im[F(k)G*(k)]sin(ujto)dk. (31) 

Using the same procedure as before (pulling out the k- independent part, 
ui(k) = luq + ui(k), and varying both the unknown angles 9 and wo*o)> one 
finds 

J fa (^Re J F(k)G*(k)e- iuJlto dk^J + flm J F(k)G* (k)e~ luJlt ° dkj . (32) 
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Again, this is equal to the standard JPD Jo - this time without dropping 
any terms in Jo due to real eigenfunctions. Therefore, the analysis from the 
previous section continues to hold for time-odd measurements followed by 
time-even measurements (and vice-versa, due to the overall time-symmetry 
of this approach). 

For more than two consecutive measurements, this problem can appear 
to be over-constrained. For example, consider a time-even measurement Qi 
at t = ti, a time-even measurement Q2 at t = and a third time-even 
measurement Q3 at t = t^. While the outcomes of Q\ and Q2 fully determine 
the solution from t\ < t < £2, the outcomes of Q2 and Q3 independently 
determine the solution from ti < t < t%. Although this implies a discontinuity 
in the solution to the KGE at t = ti-, note that it is not a discontinuity in 
4> itself, which is specified by the time-even measurement Qi- Instead, the 
discontinutity must occur at the unconstrained <f>(t — t%). 

The sort of discontinuity described in the previous paragraph is neither 
surprising nor unphysical. If <j> is supposed to describe a quantum system, 
one can only "measure" such a system by interacting with it. Because the 
interaction is outside the scope of the KGE, one can hardly hope to find 
a continuous solution without taking this interaction into account. Further- 
more, the above discontinuity is precisely where the interaction occurs, at 
t = ti. If no external system is present at the t = t<i hypersurface, then no 
information is available on this boundary, and the only external constraints 
on the system from t\ < t < t$ are the measurements at t = t\ and t = t^. In 
other words, one effectively "joins" together two adjacent regions of space- 
time into a single region - precisely as Hardy has recently argued must occur 
in any eventual theory of quantum gravity [17] . 

Incomplete measurements, however, pose a mathematical challenge for 
this two-boundary approach. A measurement that does not constrain a com- 
plete initial eigenfunction f(r) has additional "free" parameters which must 
be varied when calculating W max and W m i n . Possibly the biggest challenge 
is to determine how to implement commuting incomplete measurements, es- 
pecially when one measurement is time-even and the other is time-odd (e.g. 
L 2 and L z ). The key will be to make sure that any discontinuity at the 
second measurement does not erase the information imposed at the prior 
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(commuting) measurement. Due to the additional mathematical complexity, 
this paper can only conjecture that commuting incomplete measurements 
in this formalism will give results (nearly) equivalent to NRQM. The most 
promising path forward seems to be to ignore the time duration between 
two consecutive commuting measurements, and impose them as simultane- 
ous (but incomplete) constraints on both <fi and <f>. 

While much work is required to extend this approach to cover the full 
range of measurement theory from NRQM, such effort may be misguided 
when it comes to a possible extension to quantum gravity. That is because the 
restriction of operators to spacelike hypersurfaces will have to be revisited in 
a general relativistic framework, possibly requiring a new theory of quantum 
measurement without the use of operators at all (as in |18p. The importance 
of imposing boundary conditions on arbitrary hypersurfaces will be discussed 
in section 6. 

5 A Time-Energy Uncertainty Relation 

When examining the limits in which the previous analysis fails, one finds 
that some of these limits would manifest themselves as a time-energy uncer- 
tainty relation. Before addressing this result, it should be stressed that the 
usual position-momentum uncertainty relations are unchanged. As in stan- 
dard NRQM, a very precise measurement of <j) in position space leads to a 
large spread of <f> in k-space due to the Fourier transform ([3]). There is one 
difference; unlike standard NRQM, after a pair of complete measurements 
this approach can now reconstruct <f> between the two measurements, leaving 
no uncertainty at all (except for the free angles 9 and LUoto)- Still, this is 
true classically as well; one can "beat" Heisenberg's uncertainty principle in 
hindsight, by measuring the position very accurately at two successive times, 
and then reconstructing the intermediate velocity. So knowing a more com- 
plete solution to 4> after the fact does not contradict the usual uncertainty 
principle, which only concerns what can be known at any given time. 

In NRQM, the time-energy uncertainty principle has never been put on 
an even footing with the position- momentum uncertainty principle [19], al- 
though in a relativistic theory they must of course be intimately related. The 
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issue is not whether one can measure an arbitrarily precise energy at an arbi- 
trarily precise external clock time - in principle, one can do that in both this 
approach as well as in the standard formulation of quantum mechanics 20 . 
Here, the question becomes whether or not precise energy measurements are 
reproducible over short periods of time. 

Without using the approximation (j20|) . the exact invariant JPD (fl9|) for 
two (time-even) energy measurements can be written 



From ([33|) it is clear that if both measurements are very precise energy mea- 
surements, then they must be nearly the same energy (with the same k), or 
else the product FG* would be very small, leading to a small probability. 
But because of the delta function, there is now another constraint; the pre- 
cisely measured value of k must match up with a particular value k n (where 
k n solves uj n (k n ) = nit/to). This can only happen if to is not exactly fixed, 
but is a parameter that has some freedom to maximize W, as discussed at 
the end of section 3. This constraint on possible measurement times would 
never be directly noticeable from a single measurement because of the small 
changes in to (< 10 -20 sec) needed to set k = k n . 

But even with this freedom in the precise value of to, deviations from 
standard quantum mechanics will still occur if one cannot "ignore" the delta 
function in (|33|) . as implied by (|20p. And f|20|) is only correct in certain 
limits. The first limit derived in the appendix is (|35[) . which constrains the 
rate of change of both F(k) and G(k) for small values of k. (Specifically, for 
k < k Q: where fc = (^nm/hto) 1 ^ 2 .) For an example of such a constraint, 
consider a gaussian initial measurement F(k) = exp[—k 2 /(2a 2 )]. From (f35| . 
one finds that this corresponds to a > fc , equivalent to a spread in the initial 
energy greater than h/to- As this constraint is approached, (|33[) predicts 
that new physics will begin to emerge. The form of this new physics will be 
a failure of the ordinary probability distributions, forcing sequential energy 
measurements give slightly different results. 

Curiously, this "new physics" is almost what is expected to happen if there 
is a time-energy uncertainty relation that is analogous to position-momentum 
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uncertainty. The difference is that here the time to is not an uncertainty in 
a measurement time, but instead the duration between measurements. Still, 
in practice this might be the same thing, as it's not possible to make two 
measurements (in a well-defined order) separated by a duration smaller than 
the amount of time it takes to make the measurements. 

The other constraint derived in the appendix also concerns the rate of 
change of F(k) and G(k); see equation (|44|) . For values of k ps ko, this is 
a similar constraint to the previously discussed condition (|35[) . For larger 
values of fc, the condition becomes even less stringent, even for a precise non- 
zero energy measurement. Of course, values of k <C ko can easily violate this 
condition, but so long as the constraint (|35|) continues to hold, this will just 
be a very small perturbation on the approximation ([20)) , because the integral 
on the right side scales like k 2 , and ([35| forces the important part of the 
integral to extend out to at least k = ko. 

The conclusion is that there are many pairs of boundary conditions for 
which the approximation (|20[) does not hold good, leading to a divergence of 
J from Jo, but these pairs of boundaries also approach the time-energy uncer- 
tainty limit that is known to be experimentally difficult. Still, it is promising 
that these results yield a time-energy uncertainty relation consistent with 
known experimental limitations. 

6 Summary, Extensions, and Implications 

Although the results in this paper have not yet been generalized to the point 
where this is a full alternative interpretation of quantum mechanics, the indi- 
cations are promising that this approach might put quantum mechanics on a 
footing consistent with general relativity while also addressing many open in- 
terpretational questions. More immediately, the central result of an invariant 
scalar from which one can calculate positive definite outcome probabilities 
for single-particle solutions to the Klein-Gordon equation is itself sufficient 
motivation to further consider this novel framework. 

Three of the four postulates used in the previous sections appear com- 
patible with a general relativistic framework. These are as follows: 
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1) The correct wave equation governing generic (spinless, neutral) par- 
ticles is the Klein-Gordon equation on a complex scalar field <f>. (This is 
assumed to be a classical field, not the operator-valued field of quantum field 
theory. ) 

Although this is a postulate, it still requires some justification. As noted 
in the introduction, the KGE and SE are not equivalent in the NR limit. 
This can be seen in both the amount of necessary Cauchy data (solving the 
KGE requires twice the initial data as the SE) as well as in the doubled 
size of the KGE solution space (see Section 2). From this perspective, the 
KGE on a real scalar field would be more comparable to the SE, but such 
a field has no natural U(l) symmetry. While it is true that the solutions to 
the complex KGE in the NR limit can always be written as the sum of a 
solution to the SE and an independent solution to the complex conjugate of 
the SE, there are an infinite number of other ways to split the KGE solutions 
into the sum of two arbitrary sub-components. Assigning physical meaning 
to one of these particular sub-components under one particular splitting is 
therefore difficult to justify. For one, it raises the question of why one can 
ignore the other term; both components have positive energy, as determined 
by the zero-potential Hamiltonian H — P 2 /(2m) (as well as by the energy 
density Too of classical scalar fields). Ignoring one term also raises severe 
problems in curved spacetime where there is no covariant splitting of the 
full solution into two analogous sub-components pQ . Therefore it seems more 
natural to interpret the entire KGE solution as a single entity; that is the 
goal of this paper. 

The next two postulates are: 

2) The solution cj) is constrained by a boundary condition on a closed 
hypersurface s (with unit normal four- vector r/^), imposed by external mea- 
surements. Any infinities are dealt with by giving the mass a small imaginary 
value e, and then taking the limit e — > 0. 

3) Each possible closed boundary condition is assigned a joint probability 
distribution J = (AW) 2 , where AW is the allowed range of the real invariant 
scalar W = § s Im(4>rj^d^(j)*)ds consistent with the boundary constraints. 

The usual conditional probabilities can be generated from J via normal- 
ization. In order to compare the above postulates to a preparation-measurement 
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sequence in standard NRQM, one takes the hypersurface to be the pair of 
hyperplanes t — and £ = f o , presumably connected at some sufficiently large 
distance where the integrand in W goes to zero. Imposing boundaries on such 
hypersurfaces requires a fourth and final postulate; 

4) The measured eigenfunction of a time-even operator (from standard 
NRQM) is imposed on <f> along an instantaneous hypersurface. Time-odd 
operators are first adjusted by replacing the "i" with —ujQ 1 d/dt. 

But this use of eigenfunctions and instantaneous surfaces is the only re- 
maining tool from NRQM, and this provides a strong motivation to develop 
a measurement theory that could apply to any hypersurface. Postulate #2 
implies that any closed hypersurface will suffice, even if portions of that sur- 
face are time-like. Time-like boundary conditions are used in NRQM (e.g. 
the infinite square well), but are not imposed in the same way as measure- 
ments. This artificial distinction clearly needs to change if postulate #2 is 
correct; after all, in general relativity it is not clear whether a given surface is 
time-like or space-like until the field equations have been solved. Oeckl and 
Hardy have recently pointed this out, implying that any candidate theory 
of quantum gravity must be able to deal with regions bounded by time-like 
hypersurfaces |10j [17j . Such a generalization of measurement theory might 
also help to address some of the unanswered questions from section 4. 

The framework in this paper falls short of a full alternative interpreta- 
tion of quantum mechanics because the above analysis is only valid for the 
simplified case of a single free generic particle. Towards this end, the frame- 
work must be expanded to allow for arbitrary potentials, multiple particles, 
and non-zero values of spin and charge. While these major tasks cannot be 
accomplished in this paper, there are promising avenues for incorporating 
them into this general framework that will now be discussed. 

For non-zero potentials, NRQM simply adds a potential term to the 
Schrodinger equation. In relativistic theory, this is found to be equivalent 
to introducing an electrical potential. Unfortunately, when this type of po- 
tential is extended to the KGE, it is found that the positive- and negative- 
frequency KGE solutions respond in an opposite manner (leading to the inter- 
pretation that the negative- frequency solutions are antiparticles) . However, 
if one wants to compare the approach in this paper to standard NRQM, one 
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needs all of the KGE solutions to respond to the potential in the same man- 
ner. General relativity provides a way to add such a scalar potential: via a 
metric, <7 M „, in the weak-field limit. This does not work for NRQM, because 
the Schrodingcr equation is not covariant, but it should work for the KGE, 
which can be written g^V M V ' v <f) = m 2 4> (with covariant derivatives). Using 
this curved-space KGE to introduce potentials may be mathematically awk- 
ward, but if successful it would allow for a much easier extension to quantum 
gravity. 

For multiple particles, this approach could potentially allow a simplifi- 
cation that is not possible in NRQM: keeping the parameter space 4>(r,t) 
fixed for any number of (identical) particles. For two particles in standard 
NRQM, one expands the dimensionality of the Hilbert space, equivalent to 
a wavefunction ip(ri,r2,t). This expansion is required because of the role 
of probability in NRQM: the wavefunction ip traditionally encodes all con- 
ditional probabilities of all possible outcomes, and the number of possible 
outcomes increases dramatically with the number of particles (because dif- 
ferent types of measurements can be made on each particle). Encoding all of 
these possibilities requires a wavefunction in a large configuration space. 

But the picture of probability described in section 3 is quite different. In- 
stead of encoding all possible outcomes, the field 4> is a particular solution to 
the pair of preparation/measurements that actually happen - not all possible 
measurements. To find the outcome probability for any given experimental 
set-up, one need only compare the permitted range of W over the outcomes 
permitted by that particular set-up. Instead of the wasteful increase in the 
parameter space of ij) to deal with potential measurements that never actually 
happen, <j> only "bothers" to encode the result of the actual next measure- 
ment. This framework therefore opens up the possibility that each type of 
particle might be represented by a field in physical space, and configuration 
space would only be a tool to summarize knowledge of this physical field. 
(Such speculation agrees with the conclusion of a recent paper by Montina 
[21|.) The outcome of such a research path is far from certain, and would com- 
plicate development of a generalized measurement theory (the amplitude of 
the field must then be related to the number of particles present) but this 
also would bring quantum theory more in line with general relativity. 
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Finally, spin and charge will need to be added to this generic-particle 
framework. Incorporating spin should be straightforward, as it is known to 
involve expanding the scalar field 4> into a multi-component spinor. For exam- 
ple, the above methodology could be extended to neutral spin-1/2 particles 
simply by using the second-order van der Waerden equation instead of the 
KGE. Adding charge is less obvious, and might even seem problematic, as 
the standard interpretation of the KGE treats the b(k) terms evolving like 
exp(iujt) as antimatter. Meanwhile, the above approach requires both the 
a(k) and b(k) terms to describe a single species, consistent with standard rel- 
ativistic QM for spinless, neutral particles. But doubling the dimensionality 
of the scalar field <fi would provide parameters for both matter and antimatter 
while still retaining the above two-boundary formalism. So, for example, the 
electron/positron field would then require a four-component (Dirac) spinor; 
one doubling of the scalar field to introduce antimatter, and another dou- 
bling to permit spin-1/2. But instead of simply using the first-order Dirac 
equation, this framework requires a second-order differential equation. For 
the electron/positron field, the most likely candidate would be the square of 
the Dirac equation. In the absence of electromagnetic coupling, this simply 
reduces to the KGE on a Dirac spinor, making the connection to the above 
results quite apparent. 

This approach also has implications for quantum foundations. One on- 
going disagreement concerns whether or not the solution to the Schrodinger 
equation corresponds to some observer-indepedent "reality" , or is instead a 
construct concerning our knowledge of a system. This approach offers a mid- 
dle ground; i\) would be a construct, while the underlying system would be 
best described by some classical field (p. Without knowing the future bound- 
ary condition, the best approximation one can make concerning (j) would be 
a guess that looks something like ip, but after the fact, one can reconstruct 
what actually happened between measurements. Non-classical behavior of 
those fields, such as violations of Bell's inequality, might now be explained 
by virtue of the parameters t n , a(k), and b(k) all being hidden variables that 
are not "local" as defined by Bell [22] (because they depend on future events). 

Such a "realistic" view of <p is more plausible if e does not exactly equal 
zero, because otherwise <fi has discrete infinite coefficients at k — k n . But 
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there is some hope for a physical motivation of a non-zero epsilon; invoking 
a small but non-zero imaginary mass component has a history that goes 
back to Dirac [23] and has been revived several times since [24] [25] . Another 
promising avenue to a real <f> can be found by writing the Klein-Gordon 
equation in curved space, where the Christoffel symbols enter the equation 
in a similar manner as e |26] . Given that no quantum mechanics experiment 
has been performed in an exactly flat space, it would be possible for a slightly 
curved space-time metric to be the source of a non-zero e. 

It may be considered a disadvantage to this approach that 4> depends on 
future events. The most serious concern is the possibility of constructing a 
causal paradox. If a system <f> is constrained between two temporal bound- 
aries, then (j> contains some information about the future. At first glance, 
it seems logical that an experimenter could 1) measure <fr before the future 
boundary, 2) learn about future events, and 3) change those events to some- 
thing different. This would create an irresolvable paradox, sufficient reason 
to end this line of inquiry. But given the assumed equivalence between mea- 
surements and boundary conditions, step 1) is automatically forbidden; the 
measurement of 4> is the physical constraint corresponding to the future 
boundary condition. To measure <f> before the imposed constraint would be a 
contradiction in terms; in this framework there is no such thing as extract- 
ing information without also imposing that same information as a boundary 
condition. 

Indeed, there is a perspective that makes this retrocausal aspect an ad- 
vantage, not a disadvantage |27| . It is obvious that quantum mechanics is 
counter-intuitive, but it must be counter-intuitive for a reason - some hu- 
man intuition that fundamentally contradicts some physical principle. One 
example of this would be the well-known conflict between our directed ex- 
perience of time and the more symmetric treatment of time in fundamental 
physics. If the counter-intuitive aspects of quantum mechanics could be ex- 
plained via classical fields symmetrically constrained by both past and future 
events, then it would be a mistake to reject such a solution based solely on 
our time-asymmetric intuitions. 

Regardless, the ultimate test for any theoretical proposal is experiment. 
Perhaps the strongest reason to pursue this line of research is that it leads 
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to results which contradict standard quantum mechanics (albeit only in 
experimentally-difficult regimes). Hopefully an incorporation of potentials 
and a relativistically correct measurement theory will lead to practical ex- 
periments which can confirm or deny the validity of this general approach. 
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This appendix derives the limits in which the approximation (|20)l holds good. 
Recall k n are the solutions to the equation u>„(k n ) = nir/to, and no is the 
smallest value of n for which there is a solution. It will also be useful to define 
the quantity fc„+i/2 using the similar equation Lo n {k n +i/2) — (n + l/2)7r/£o- 
To compare to NRQM, one is interested in the limit where u> « mc 2 /% + 
hk 2 /(2m). 

Working backwards from the right side of (f2T7|) . using spherical k-coordinates, 
the "radial" component of this integral can then be written as the discrete 
sum of integrals 

/>oo 

/ Re[FG*]cos{wtu)k 2 dk = ££L„ +i f^j*Re[FG*}cos{ujto)k 2 dk + 

j Q kl/2 Re[FG*}cos(ujto)k 2 dk. (34) 

The final integral (where fc„ 0+1 / 2 is written as k\/2 for clarity) is different 
from the others because the lower limit must be forced to zero, and the precise 
value of k no _i/2 depends on the values of m and to - indeed, it may not even 
exist. In the special case that fc„ -i/2 — 0, then fc„ 0+1 / 2 — (inm/hto) 1 ^ 2 = 
ko- This value, ko, is a reasonable estimate for the value of k no , given that 
m and to are not known precisely. 

The last integral in (|34|) is just one of many similar integrals that scale like 
k 2 , so this lowest-k integral it is not likely to be important unless F(k)G*(k) 
is only large for values of k < fco- But such a narrow range is only possible 
if either (1/ F)dF / dk or (l/G)dG/dk evaluated at k < ko is on the order of 
ko- This leads to the first constraints on the approximation (|20p . 
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This constraint is discussed above in section 5 on the uncertainty principle. 
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Each integrals in the sum in (f34|) will be referred to as Y(k n ); writing out 
the It-dependence explicitly (in the NR limit) 



Y(k n )= I RelFi^G'^k^cosiuJoto^osi^k^- 



sin^uJo^sin^^^^dk, (36) 



using the constants ljq = mc 2 /% and a — y / ht Q /(2m). 

Taylor expanding the nth integrand around g(k„) = Re[F(k n )G*(k n )], 
the zero-order term in the expansion can be integrated exactly: 



Y(k n ) = g(k n )—sin(uJoto + a 2 k 2 ) + 0[g'(k n )]- 



k 

2a 2 ' 

l( k ri)\ -3— (5 [cos(uj t )S(ak) + sin(u to)C(ak)} 



■ (37) 

Here C(x) and S(x) are the Fresnel integrals, which for x > ir can be excel- 
lently approximated by [28j 

C{x) = — H — ; sin(x 2 ) == — cos(x 2 ) + 0(x~ 5 )sin(x 2 ) + 0(x~ 7 )cos(x 2 ), 

(38) 

S( x ) = 7: -^^cos(x 2 ) sin(x 2 )-0(x~ 5 )cos(x 2 ) + 0(x~ J )sin(x 2 ). 

2 \/2lTX V87TX 3 

(39) 

With x — ak, the constraint x > tt corresponds to k > kg, reinforcing the 
earlier discovery that this approximation fails when the first integral in (|34|) 
dominates the others, and therefore is equivalent to the earlier constraint 
([33)) . Plugging in the limits in (f3"T|) , dramatic simplifications occur because 
sin(uj n+1/2 t Q ) = (-1)" and cos{uj n+1 / 2 t a ) = 0; 



Y(k n )= g(k n )(-iy 



(fcn+l/2 + fcn-l/2) ( fc n+l/2 + ^n-l/2) 
2^2 + 8^ 

O(a- 10 A:- 7 ) + O(g'). (40) 

r(fc„) « g( kn )(-ir^ + O(a- 10 k~ 7 ) + O(g'). (41) 
Mo 

The last step is an approximation, good to better than 1% even for the first 
term (n = no + 1), and rapidly improving for the higher terms in the sum. 
Within a constant, the first term in (|4"Tj) is exactly the left side of ([20]) . with 
the correct scaling of both k n and to- 
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The next term in the Taylor expansion, g'(k) — d Re[F(k)G* (k))/dk, is 
unimportant if the variation with k is sufficiently slow. This integral is not 
much harder than the last, because the k 3 cos(k 2 ) integral can be done by 
parts without requiring the Frcsncl integrals. To highest surviving orders, 
this term becomes 



0(g')=-g'(k n ) 



( 1 )" I 1.2 , L.2 

/2 + K n-l/2 



2a 2 
k. 



-h 2 

\ K n+l/ 

kn+1/2 + fcn-1/2 + iKll/2 + fc r7-l/ 2 )/( 4 « 4 )] } . ( 42 ) 



Y{k n )^g(k n )(-l) 



n 2mk n it 2 - 8 
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kl 



-O(a- 10 k- 7 ). (43) 



Again, this last step is a quite accurate numerical approximation that 
increases in accuracy with n. Comparing the magnitude of these two terms, 
one finds that the approximation (j20|) fails unless the fourier expansion of 
both boundary conditions vary sufficently slowly over k such that 



^[F{k)G*{k)\ 



«170|F(fc n )G*(fc„)|§ 



(44) 



This result is also is discussed above in section 5 on the uncertainty principle. 
The conditions on the higher order terms in the Taylor expansion are not 
worked out here, but they should also be taken into account. 
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